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In the interaction of high-power laser beams with solid density plasma there are a number of 
mechanisms that generate strong magnetic fields. Such fields subsequently inhibit or redirect electron 
flows, but can themselves be advected by heat fluxes, resulting in complex interplay between thermal 
transport and magnetic fields. We show that for heating by multiple laser spots reconnection of 
magnetic field lines can occur, mediated by these heat fluxes, using a fully implicit 2D Vlasov- 
Fokker-Planck code. Under such conditions, the reconnection rate is dictated by heat flows rather 
than Alfvenic flows. We find that this mechanism is only relevant in a high /? plasma. However, the 
Hall parameter w c T e i can be large so that thermal transport is strongly modified by these magnetic 
fields, which can impact longer time scale temperature homogeneity and ion dynamics in the system. 


Understanding the O (lO 2 ) T magnetic fields that 
can develop in high-power-laser interactions with solid- 
density plasma m is important because such fields sig¬ 
nificantly modify both the magnitude and direction of 
electron heat fluxes |B]. The dynamics of such fields evi¬ 
dently has consequences for inertial fusion energy appli¬ 
cations cm, as the coupling of the laser beams with the 
walls or pellet and the development of temperature in¬ 
homogeneities are critical to the uniformity of the implo¬ 
sion. There is a significant interplay between heat fluxes 
and magnetic fields: in semi-collisional plasmas heat flux 
can be the dominant mechanism for transporting mag¬ 
netic fields in addition to currents or bulk ion flow JTO] . 
This effect, arising due to an electric field analogous to 
the Nernst-Ettingshausen effect in metals El, has been 
shown to be significant in laser heated plasma C2HI3. 
The Nernst effect in plasma uni arises as a consequence 
of the velocity dependent collision frequency of electrons 
in plasma. Since the faster, “hot”, population of elec¬ 
trons are essentially collisionless, the magnetic field is 
“frozen” to them, whereas the collisional, “cold”, por¬ 
tion of the distribution function is able to diffuse across 
field lines. Hence, magnetic fields can be advected with 
close to zero net current by “hot” electrons. 

In heating plasma with a finite laser spot, an azimuthal 
magnetic field about the heated region arises through 
the Biermann battery effect m For multiple spots in 
close proximity, as in inertial fusion, these magnetic fields 
will be in a configuration with oppositely directed field 
lines. Under such conditions, magnetic reconnection of 
field lines may be expected to arise. Magnetic recon¬ 
nection has been intensely studied in space plasmas, but 
more recently, laser inertial fusion relevant scenarios have 
been investigated [DBMS] . 

In Sweet-Parker theory, a resistive region between 
plasma inflows with resistivity 77 allows magnetic field 


lines to diffuse and change topology, leading to jets out¬ 
flowing at Alfvenic speeds, va El [201 . However, ob¬ 
served reconnection rates are rarely well described by this 
model. In case of a small diffusion region L for which the 
Sweet-Parker width, Ssp = L/^JS where S = vaLro/ii 
is the Lunquist number, is smaller than the ion inertial 
length, c/ujpi , Hall physics is relevant and reconnection 
is no longer resistivity dominated. Rather, augmented 
by the inclusion of Hall physics, reconnection rates have 
been shown to be significantly faster suggesting that the 
dynamics at such small scale lengths contribute strongly 
to reconnection. Recently, researchers have been inter¬ 
ested in the intermediate regime between collisionless and 
collisional reconnection [2fH24j . where the Nernst effect 
is important. Given the analogous form of the Nernst 
term and Hall current term in Ohm’s law, it is natural to 
assume that the Nernst effect may enable reconnection 
in a similar manner to Hall reconnection, but with the 
electron currents replaced by heat fluxes. 

In this Letter, we show that under conditions similar 
to those found in hohlraums, where heat flux effects in 
Ohm’s law are important E, reconnection of field lines 
can occur. The heat fluxes that are generated by the laser 
hot-spots drive reconnection through advection at the 
“Nernst” velocity vp- The Nernst effect allows magnetic 
field advection without an associated electron current, 
which is different than the standard Hall effect within the 
reconnection layer; this breaks the Alfvenic constraint (at 
least within the parameters considered) and allows char¬ 
acteristic reconnection rates of E z /(B 0 vt) rather than 
E z /(Bqva)- We show that this can occur for conditions 
described by a dimensionless number describing the ra¬ 
tio of Nernst to electron flow velocities. We find that 
this mechanism is only relevant in a high /3 plasma, i.e. 
where the ratio of thermal pressure to magnetic pressure 
is large. However, the Hall parameter w c r e j can simulta- 
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neously be large so that thermal transport is strongly 
modified by magnetic fields, which can impact longer 
time scale temperature homogeneity and ion dynamics. 

The Vlasov-Fokker-Planck (VFP) equation, is solved 
together with Ampere’s and Faraday’s Laws to describe 
the plasma. The code we use, Impacta 25,, 26], uses 
a Cartesian tensor expansion gzi, with the distribution 
function expanded as f(t, r, v) = /o + fi-v + / : vv + ... 
This expansion can be truncated in a collisional plasma, 
as collisions tend to smooth out angular variations in 
the distribution function, resulting in a close to isotropic 
distribution in the centre of mass frame, represented by 
fo. Higher orders are successively smaller perturbations, 
/ 2 <C fi <C fo etc. Using the Lorentz gas approximations, 
electron-ion collisions appear in the equation describing 
the evolution of fi in the ion center-of-monrentum frame 
as an effective collision frequency oc 1/v 3 : 


dfi 

dt 


+ vV/o — 


_ 2 _ fL 

5v 3 dv 


eE dfo 
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where Y = 4Tr(e 2 /ATTeom e ) 2 lnA e i. In Impacta, terms 
up to and including / are retained. 

In reconnection studies, Ohm’s law is of crucial sig¬ 
nificance. We can formulate a generalized form of 
Ohm’s law for this velocity dependent collision opera¬ 
tor by multiplying by v 3 and taking the current moment 
( 47r / 3 ) /o°° ■ ■ • v3dv - 
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where the effective resistivity is 
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the magnetic convection velocity by heat flow m is 


v T 


(v^ 3 ) , j 

2(v 3 } en e 


the inertial term (d/dt) is neglected, valid for a suffi¬ 
ciently collisional system, a term contracting with E, 
(vvu)/(2(v 3 )), is assumed to be small, and velocity mo¬ 
ments are defined by 

(v n ) = £ fo°° fo v n+2 dv , 

(w n ) = it fo°° f i v n+3dv - 
(w^) = jStr j2V n+4 dv . 

The last two terms in Eqn. [2] combined play the role of 
the pressure tensor term normally used in Ohm’s law. 


To express Eqn. [2] in a more familiar form, if we assume 
that the distribution function is a Maxwellian speed dis¬ 
tribution multiplied by a function of angle only, it can be 
shown that Eqn. [2] reduces to 


E = m + 
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where P_ is the full electron pressure tensor and we have 
neglected a ~ (vv)/ (v 2 ) correction to the VT e term. We 
do not use Eqn. [3] in this study, but instead compare 
the results from the Vlasov-Fokker-Planck code with the 
more general Eqn. [2] 

To gain some insight into the physical meaning of vj-, 
Haines showed, using a model 1/v 2 collision operator [TP] , 
that the Nernst velocity could be related directly to the 
heat flux by v-r ~ 2q e /(5p e ) and how it relates to terms 
in Braginskii’s equations [6]. 

To parameterize under what conditions the situation 
we describe may occur, we can compare the relative mag¬ 
nitudes of the Hall term, j x B/en e , and the heat flow 
term, vy x B in Ohm’s law to generate a new dimension¬ 
less number: 


H n = 


en e v T 


1 K c , 


5 U> r T e , 


) = \k C l f3u c T e i (4) 


where co c T e i is the Hall parameter, k[ j_ is the normalized 
perpendicular thermal conductivity coefficient p5 : and 
we have used the heat flux component q eJ _ ~ n±'VT e to 
estimate v T and assumed the gradient scale lengths for 
the temperature and magnetic field are similar. The nor¬ 
malized skin depth, S c = c/(vthW pe r e ,:) serves as an inde¬ 
pendent parameter in Eqn. [4] A small skin depth relative 
to the mean-free-path means that electron currents are 
inhibited, but the senri-collisional behavior still allows for 
electron energy transport. /3 is the ratio of thermal pres¬ 
sure to magnetic pressure. In the limit of large w c r e j, the 
k c _ j_ approaches the asymptotic limit = 7 [/(uj c T e i) 2 , 
where 7 ^ is a coefficient between 3.25 and 4.66 depending 
on Z PB], Hence, for large w c T e j, = ( 7^/5 )/3/uj c T e i 
and can therefore only be significant for a high ft plasma. 

An important parameter in magnetic reconnection is 
the Lundquist number S. We can also introduce an 
analogously formulated Nernst-Lundquist number, Sn = 
VTLyLo/r], which is defined according to the usual defi¬ 
nition, but replacing the Alfven velocity with the more 
relevant Nernst velocity. The relationship between these 
two dimensionless parameters is 


S]\J — S \r(jJcTei 



( 5 ) 


From these dimensionless numbers, we can see that for 
an interesting heat-flux reconnection problem (i.e. for 
u) c T e i > 1) dominated by Nernst effects ( Hn 1), the 
Nernst-Lundquist number must also be large, Sn 1. 
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This means that resistive effects will be small, and there¬ 
fore anisotropic pressure-like (/ 2 ) effects must be in¬ 
cluded in Ohm’s law to support the electric field at the 
X-point. In ref. [ 22 ] , Daughton et al included heat flux 
effects in their reconnection study, but for their system 
Hn ^ 1 so the thermal contribution was small. Here 
we examine a situation where Hn 1 where heat flux 
effects dominate. 

As in ref. [251 [26], we use a normalization scheme with 
time normalized to r n = kr e il?>\pK and velocity normal¬ 
ized to t>tho = v^/caTeo/we- As a result, lengths are 
normalized to the electron mean-free-path A m f p = VthX n - 
The computation is performed in a domain defined over 
the range — 100A m f p < y < 100A m f p and — 1500A m f p < 
x < 1500A m f p . The cells near the boundary in x is expo¬ 
nentially increasing in step size such that they can be con¬ 
sidered “far away”. The domain of interest in x, where 
the cell size is constant, is —400A m f p < x < 400A m f p . 
The numerical resolution in the runs shown in the paper 
is Ax = 13.3333A m / p , Ay = 3.125A m fp, Av = 0.0625 v t h- 

The connection between the normalized quantities and 
real parameters is made through the ratios Wth/c and 
u} pe T n . Here, Uth/c = 0.08 and Lo pe T n = 125 are chosen in 
order to put the system into inertial confinement relevant 
conditions, corresponding to a temperature T e0 = 1.6 keV 
and electron number density n e = 2.5 x 10 22 cm -3 . A 
magnetic field of B 0 = 1 corresponds to a field strength 
of 400 T (4 MG). 

The magnetic field is generated through the Vn e x 
VT e mechanism. We introduce an out of plane 
plasma density gradient of the form dn(x,y)/dz = 

(j—( x / r o) ((y-tymax)Ao)~ _j_ g — {{y— i/maxlAo)"^ where 

L n = 50 and tq = 50, by adding a 2 component of electric 
field. This gradient is switched off at t = 800r„ to prevent 
excessive magnetic field generation. The temperature 
profile is accomplished by heating the plasma near the y- 
boundaries of the system using an inverse bremsstrahlung 
heating operator [29] with a profile H(x,y) = 
HqC~ (* c / r o) ((y-tymax)/ro) ((y—ymax)/ r o)~'j where 

Ho = 0-5, corresponding to a laser of intensity 2.5 x 10 14 
Wcm -2 . The heated regions result in strong heat fluxes 
in the y direction that advect the magnetic field lines 
inwardly towards the reconnection region. Ions are sta¬ 
tionary in the simulation to isolate these effects, which 
may be justified physically in the case of the walls of a 
hohlraum as they are heavy ions (for example gold) [30]. 
Simulations run with ion motion show similar behavior. 
A thorough study of the Nernst and bulk flow advection 
of magnetic fields is in Ref. [Ill- 

Figure [T] shows output from the simulation at a time 
19000 r n . (a) shows the magnetization of the plasma, B 
and (b) illustrates the temperature profile of the system. 
The Nernst velocity (c) is approximately 10 2 larger in 
magnitude than the maximum current (not shown). The 
flow direction of the Nernst velocity (calculated directly 



x / A m f p 


FIG. 1: At a time t = 19000 r„ into the simulation (a) 
B/m^ ei /e, (b) T e /T e0l (c)v T /u th0 in the x - y plane, 
(d) B /mv e i/e and (e) T e /T e0 at t = 27000 r„ 
Note: The axes are not square 


from the distribution function) indicates that thermal en¬ 
ergy is being brought inwards in the y-direction towards 
the reconnection region and is subsequently redirected 
outwards in the x-direction, carrying the magnetic field 
with it. Distinct “jets” of heat flux are formed out of the 
reconnection region. 
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Figure [TJd) shows the magnetic field profile after the 
majority of the flux has reconnected, at a time 27000 r n 
into the simulation, which corresponds to approximately 
0.6 ns. The reconnected field lines are then advected by 
the Nernst jets towards the x boundaries. Figure [lje) 
shows the temperature profile at the same time as the 
magnetic field in (d). The outward heat flow in x from 
the reconnection process causes the change in the tem¬ 
perature profile from (b) to (e). 



y / A m fp 


FIG. 2: Illustration of the contribution of the different 
components of Ohm’s Law in Equation ([2]) taken from 
the simulation at a time t = 11000 r„. 

(a) E z calculated from the code 

, (b) fjj z , (c) [j x B] z , (d) [v T x B]-,(e): 

(f) Sum of all contributions (b-e) 


V-(vvu 3 ) 

2(i>3) 


the magnetic field near the reconnection region. Fox et. 
al. [23. account for this effect in the calculation of the 
local magnetic field, and we perform the same correction. 
We find that in our simulation, E z /Bvt ~ 0.1, as shown 
in fig.[|[a). 




t (r n ) 


FIG. 3: (a) The reconnection rate (b) /? - the ratio of 
thermal pressure to magnetic pressure (c) w c r e j as 
functions of time. 


The quantity E z is the rate at which magnetic flux 
crosses the neutral point. In the case of oppositely di¬ 
rected magnetic fields, B x , the reconnecting magnetic 
field, By , is generated through Faraday’s Law by the out 
of plane electric field, dE z /dx in a 2-D Cartesian geom¬ 
etry. We can analyze the various contributions from the 
generalized Ohm’s Law, eq. pi) by directly calculating 
the velocity moments. Figurcpl shows the out-of-plane 
electric field, E z , and four of the terms that contribute 
to it. Anisotropic pressure tensor-like terms almost en¬ 
tirely support E z at the X-point where the flows diverge, 
with a small contribution from the resistive term. The 
vj x B term provides an analogue of the Hall current, 
with the actual Hall current j x B being negligible. The 
sum of just these moments of the numerical distribution 
function agrees well with the electric field taken from the 
code (which in these calculations includes electron iner¬ 
tia). Using the terms in eq. ([3]) instead, similar results 
are obtained, with the small difference being due to the 
non-Maxwellian distribution that develops in the recon¬ 
nection region. 

By convention, as in I23i . the reconnection rate co¬ 
efficient is reported as E z /Bva where va is the Alfven 
velocity and typical rates associated with fast reconnec¬ 
tion are E z /Bva = 0.1 ~ 0.2. In our simulation, the ions 
are fixed, and consequently, Alfvenic flows are nonex¬ 
istent. The characteristic flow velocity for the flux is 
clearly vj-. There is a marked increase in the strength of 


Figure |3])b) illustrates the evolution of /3, the ratio of 
thermal pressure to magnetic pressure. The sharp peak 
arises due to rapid heating of the plasma, and then the 
subsequent decrease comes from the compression of the 
magnetic field flux before the reconnection process can 
begin. Once the anisotropic pressure-like term, fig. |2](e) 
supports the out-of-plane electric field, E z , across the re¬ 
connection layer, the field compression is maintained and 
eventually reduced, while the plasma is heated due to the 
decrease in transport inhibition because of the reconnec¬ 
tion process. This corresponds to the steady increase in 
/3 as observed after t = 13000r n . Figure |3](c) illustrates 
the magnetization of the plasma over time. The initial 
rise in uj c T e i is due to the compression phase of the mag¬ 
netic field. After this period, a plateau arises because 
while the magnetic field decompresses due to reconnec¬ 
tion, the plasma heats in the reconnection region, effec¬ 
tively increasing r e i. The steady decrease in the late-time 
behavior is attributed to the magnetic field decompres¬ 
sion as the reconnected field lines relax from the Nernst 
outflows. 

When simulations are performed with different values 
of uip/uei and Vth /c, i.e. different plasma density and 
temperature, they evolve similarly for fixed 5 C , with a 
reconnection rate stabilizing close to E z /Bvt ~ 0.1, as 
expected from Hn, which only depends on S c for fixed 
uj c T e i. Hence, the ratio of the skin depth to the collision 
mean-free-path is the important consideration for this 
mechanism. One significant difference with Hall recon- 
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nection that we wish to highlight is that quasi-neutrality 
can be maintained throughout the system and therefore 
there is no necessity for ion motion outside of the recon¬ 
nection region to maintain dynamic equilibrium. Redi¬ 
rected heat flows by magnetic reconnection can result in 
a redistribution of thermal energy and reconnection of 
field lines can remove thermal transport barriers. Since 
strong heat flows and magnetic fields are expected in the 
interior of hohlraums, understanding this mechanism can 
be expected to be important for inertial fusion energy, in 
particular because reconnection may mitigate the ther¬ 
mal transport inhibition by magnetic fields that could 
affect the uniformity of the drive. However, the mag¬ 
netic reconnection could also lead to the production of 
energetic electrons. 
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